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We use random matrix theory (RMT) to study the first two moments of the wave power transmit- 
ted in time reversal invariant systems having ergodic motion. Dissipation is modeled by a number 
of loss channels of variable coupling strength. To make a connection with ultrasonic experiments 
on ergodic elastodynamic billiards, the channels injecting and collecting the waves are assumed to 
be negligibly coupled to the medium, and to contribute essentially no dissipation. Within the RMT 
model we calculate the quantities of interest exactly, employing the supersymmetry technique. This 
approach is found to be more accurate than another method based on simplifying naive assumptions 
for the statistics of the eigenfrequencies and the eigenfunctions. The results of the supersymmetric 
method are confirmed by Monte Carlo numerical simulation and are used to reveal a possible source 
of the disagreement between the predictions of the naive theory and ultrasonic measurements. 



I. INTRODUCTION 

The statistics of waves in complex disordered and ray- 
chaotic structures have been well modeled in recent years 
by random matrix theory (RMT). The eigenstatistics of 
such structures are ergodically equivalent to those of cer- 
tain classes of random matrices. This has been estab- 
lished by an enormous literature, both experimental and 
theoretical, and appHes to the complex dynamics of com- 
pound nuclei and also to the somewhat simpler case 
of closed nondissipative wave billiards with chaotic ray 
trajectories 0- More recently attention has turned to 
the case of structures with open loss channels and/or in- 
ternal dissipation [3-8]. This evolution of focus has been 
dictated by the physical systems available, for which it 
is difficult or impossible to ehminate absorption, and/or 
minimize the dissipative effect of the channels used to 
launch or detect the waves. 

Of the many statistics that could be considered for 
such structures, perhaps the simplest experimentally ac- 
cessible one is the relative variance of the power trans- 
mission. This quantity is related to cross section fluctu- 
ations in nuclear reactions, it is accessible in microwave 
experiments, and is of long standing interest in acous- 
tics |M ■ In R-ef ■ (§1 relative variances measured in a 
dissipative 3D ultrasonic billiard were compared to the 
predictions of a simple theory which assumes, that the 
eigenstatistics are identical to those of the nondissipative 
case. Such an assumption is strictly true only for very 
special cases of dissipation. The theory was found to con- 
sistently overestimate the relative variance of the mean 
square transmitted signal. Our chief interest here is to 
develop a more rigorous theory for that variance, and to 
compare it with the predictions of the naive theory. 

As an adequate theoretical model of such a structure 
we will use an effective random matrix theory description, 
with a random matrix H replacing the wave equation's 
linear differential operator (see and references therein 



for a more detailed discussion). Then the complex am- 
plitude of the transmitted wave is simply proportional to 
the ijth matrix element of the resolvent associated with 
the wave equation: G{E) = [E I + ie I — H + ir]~ . 
Here, / is the identity matrix, the matrix F describes 
coupling to internal local- in-space dissipative channels, 
the parameter e > describes uniform losses, and E 
is the spectral variable. The same model describes mi- 
crowave billiards, ultrasonic billiards and nuclear reac- 
tions. The real symmetric positive semi-definite loss op- 
erator r can be written in terms of the states of the chan- 
nels (r = W^W in definitions of Ref. illj) or in terms of 
absorption mechanisms. It is generally taken to be only 
weakly dependent on E. Thus both open and closed dis- 
sipative systems are described by the same model. When 
losses are negligible the experimental systems are usually 
invariant under time reversal. The appropriate choice for 
the corresponding random matrix H should, therefore, 
be a real symmetric matrix taken from the Gaussian Or- 
thogonal Ensemble (GOE). 

Our quantity of prime interest is T = jGy (E) \^ ,i j, 
i.e. the product of retarded and advanced Greens func- 
tions (propagators) : Gfi{E) = [E I + ie I - H + iT]~.^ 

and = (Gj^(i?)) respectively. Except for slowly 
varying factors of receiver gain and source strength, the 
quantity T represents the ultrasonic power of Ref. 

For general non-perturbative statistical studies the 
only generally applicable tool known at present is re- 
duction to the Efetov's zero dimensional supersymmet- 
ric non-linear ct- model In this way the problem of 
calculating RMT ensemble averages reduces to perform- 
ing a definite finite-dimensional integral over a space of 
supermatrices. The zero dimensional cr-model can also 
be derived from the assumptions of RMT and is 
applicable to a variety of quantum-scattering problems 
formulated in terms of random Hamiltonians [8-12]. Ide- 
ally, once the quantity of interest is expressed in terms 
of the products of resolvents of the effective Hamiltonian 
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H — ir, its mean, its variance, and sometimes its distri- 
bution function can be obtained. 

Technical details of the supersymmetric reduction pro- 
cedure depend essentially on the basic symmetries of the 
underlying ensemble. It is well known, that working with 
the orthogonal ensemble leads to calculations, which are 
more technically involved than those of similar calcula- 
tions for systems with broken time-reversal invariance 
(TRI). In the latter case, the corresponding ensemble is 
the Gaussian Unitary (GUE), and one can go as far as 
calculating the full distribution function of transmitted 
wave power in ergodic systems Q. Unfortunately, the 
existing experimental results on power transmission are 
only available for systems with preserved TRI. 

The main goal of the present work is to explore trans- 
mitted power statistics for ergodic TRI systems. We find 
that it is possible to derive explicit analytical expressions 
for the first two moments of this quantity. We wish in 
particular to explain the differences, seen in between 
the predictions of the oversimplified ("naive" perturba- 
tive) theory and experimental measurements. 

In section II we use the supersymmetry method to de- 
rive expressions for mean and variance of transmitted 
power. In Section III, we confirm the results by numer- 
ical Monte Carlo calculation, and compare them with 
the results of the perturbative method of Refs. 0, 
In Section IV we investigate a hypothesis to explain the 
longstanding discrepancy between lab measurements in 
ergodic acoustic systems and naive RMT predictions. A 
summary is given in Section V. 



II. SUPERSYMMETRIC CALCULATION FOR 
MEAN AND MEAN SQUARED POWER 

A. The system 

In an ergodic system characterized by a random NxN 
Hamiltonian H and a dissipation matrix F, a matrix el- 
ement of G {E) = {EI + iel — H + iV)^^ represents the 
response spectrum (with E being the spectral variable). 
Its squared absolute value Gij [E) G*j (E) denotes the 
spectral power density. 

The elements of the random matrix H are zero- 
centered Gaussian variables, and because we deal with 
with power transmission inside time-reversal invariant 
systems, the matrix H is real and symmetric. The rele- 
vant random matrix ensemble is, therefore, the Gaussian 
orthogonal ensemble (GOE). Because of the inherent or- 
thogonal invariance the dissipation matrix may be chosen 
to be diagonal: T — diag{7,7, ...7,0, ...0}, as we always 
can express our matrices in F's natural basis. The num- 
ber Af < iV of nonzero entries can be interpreted either 
as a number of equivalent open channels in the scattering 
system or, as a number of equivalent localized "dampers" 
in a closed system with losses 0. Note, that conver- 
gence generating parameter e > can be interpreted as 
the coupling to infinite number of external channels or. 



as uniform dissipation. 

We are interested in the statistics of the wave power 
T transmitted from a source at point j to a receiver at a 
different point i inside the system: T = Gij (E) G*j (E) , 
1 ^ j ^ N j ^ i (no summation over i and j) 0. 



B. Basic definitions and identities 

To obtain expressions for the first two moments of the 
transmitted power T, (the bar indicates the ensemble 
averaging) we adopt a procedure similar to that of Ref. 

We start with constructing a generating function Z 
for our quantities of interest by introducing 4-component 
supervectors ^p-. 

$J = {xf , Xj, Spilf, Spi2f} , p - 1, 2 , 

where the components of TV-dimensional vectors S are 
real commuting variables, the elements of the vectors 
X are anticommuting variables (Grassmannian), and 
stands for the vector transposition. Index p is used 
to distinguish between retarded {p — 1) and advanced 
(p = 2) Green's functions. The latter can be ob- 
tained from the generating functions: Zp{E,^p) = 
/ [d$p] exp{(i/2) £p (i?, $p,Olp)}, where the 'actions' £p 
are defined as: £p {E, ^p,Zp) = (S)p + 3p) ^p in terms 
of the block-diagonal 4x4 symmetric supermatrices [7- 
10]: 

Dp ^ {EI ~ H)® Lp + i {el + T)® KpLp, 
Li = diag{/2,/2} ,^2 = diag {12,-/2} , 
Ai = diag{/2,/2} , A2 = diag{-/2, -I2} , 

ai = diag {0,0, J(i) , J(2) I , 32 = diag {0, 0, J(3) ^ j(4) I ^ 

Here j'^^^ are NxN symmetric source matrices, and 
the integration measure is just a product of indepen- 
dent differentials of commuting and anticommuting vari- 
ables. The generating function for the power moments 

T ^ Sy'ST/ and = (D-.^D*-/)" then can be 
shown to have the following representation: 

z{E,d) (i;,ai) ^2(^,^2) 

= J [d$]cxp|^£ (i?,$,3)|, (1) 

in terms of 8 x 8 block-diagonal supermatrices 
D = diag {Di, 2)2}, L = diag{Li,L2}, A = 
diag{Ai, A2}, a = diag 32} and $ = {«'i,$2}, 

£(/;,$,a) = £i(^,$i,ai) + £2 (s, $2,^2) = 
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The Gaussian integral over the supervectors in Eq. (QJ 
can be also written as a superdeterminant 

ZiE,Z)= n ZpiE,^p)= Yl Sdet -1 (Dp + ap) . 

Differentiating this expression with respect to elements 
of the symmetric source matrix 3 one finds (cf. [Tllll^l: 



(2) 



i^2[<I>] = ^ (1)1,5(1)2.^(2)1, 5 (2)2, 



Now, we proceed with GOE averaging of the above ex- 
pressions for the moments of the transmitted power. In 
what follows we use the overbar to denote the averag- 
ing over H with the weight exp|— (^N/Av^)TrH'^H^, 

so that H^jHki = {v'^ /N) {S^kSii + SiiSjk), i-e. the en- 
semble averaging. It can be performed exactly with the 
help of the identity: 
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(3) 



relating both T and to the Gaussian integrals over the 
supervector components. Using the shorthand notation 
(...)^ = / (. . .) exp{i£ {E, $) /2}, we can write 



T = (Fi [$]) 



(4) 



exp <^ -$t (i/ (g) L) $ I ^ exp 



-StrA^ 



AN 



where we introduced a new 8x8 supermatrix: A = 
Eili The elements of the supermatrix A 

are labeled as follows: 



= {F2 [$]) 



(5) 



A 
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where we introduced the following products of the com- 
muting components of the supervectors 



where m,n — 1, . . . , 4.. With the help of these notations 
we can express T and in a unified form via the repre- 
sentations 



M$] Fi^2 [«>] exp { ^-E<^^L-^> - (r ® A) L$ - ^Str^^ - |strAA 



as both formulas differ only in the form of preexponent 
factors F. 



C. Performing ^-integration 

The next step of the supersymmetric calculation is the 
so-called Hubbard-Stratonovich decoupling (ifl ll^ : 

exp l-^StrA^ - IstrAAl 

[dR] exp <^ - — Stri?2 + i—NSivRK + i-^ivRA 



{Fi,2 J [dR] exp |-^Stri?2 + i^NStiRA 

X J [d^]F,^2mexp!^^'-<i>^L'/^r'L'/HY (6) 



where we defined 8A^ x 8A^ supermatrix /: 

/ = [-EI (g) h - vl(g) R - i {T (g> A)]~^ 
= [{In ® 4 - «r ® (A0-1)) ©] , 

with = —Els ^ vR. In Eq. Q we can integrate out 
<i>— variables, using Wick's theorem for supervectors, and 
bring the remaining integral into a form suitable for a 
saddle point approximation in the limit N 00 [ifiLIT^ . 
Then for i,j > M we obtain: 



}• / 



[d$]fi,2 [$]exp|-^$t;-i$| 

= Fi,2 (Sdet/)^/' . (7) 



Here we introduced the notations 



Fi[^-']=\[i^-X + i^-')l}'^ (8) 
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and 
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where 6^^ = jc^^ + /2. At this point we 

summarize the results for T and separately: 

T = J [dR] Fi [e-^] exp {~NC [R] + S£} , (10) 
T2 = y" [dR] F2 [25-^] exp {-NC [R] + S£} , (11) 



where the exponential is given by 

C [R] = istri?2 + istrln {-Eh - vR) , 



(12) 



SC ^i—NStrRA 
2v 



M 



Strln 



Is - ijA {~EIs - vR) ^ 



(13) 



The remaining step is to carry out integration in Eqs. 
Hl()lll|l by the saddle point method in the limit of large 
N. The stationarity condition for C [R] yields the sad- 
dle point equation Rs = v/ {—Els — vRg). Its solution 
is given by a saddle-point manifold in a space of 8 x 8 
supermatrices (icL lll]|: 



Rs = 



E 
2^' 



E 
2^ 



h - TTWQ. (14) 



Here u denotes the mean eigenvalue density given 
f or the GO E by the Wigner semicircular law v — 
\J ^v^ — E'^ I (27ru^). After integrating out the massive 
Gaussian fluctuations around the saddle point manifold 
in Eqs. lllOlllll . the first two moments of the transmitted 
power are expressed as integrals over the supermatrices 



T-iAT 



{-KV)' 



[pressed 



;Sdet-^^/2 



X exp <j --eTTi^iVStrQA 
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= (^2 [Q])q . 



(15) 



(16) 



This step completes derivation of the zero-dimensional 
nonlinear cr-model. 

D. Performing Q-integration 



To evaluate the superintegrals in Eqs. Ijl5ll6|l . we need 
to calculate Fi [Q] and F2 [Q] first. At this point we em- 
ploy the Verbaarschot-Weidenmiiller-Zirnbauer (VWZ) 
parameterization for the matrix Q. Both Fi [Q] and 
F2 [Q] are the functions of matrix elements of Q, obtained 
by the formal substitution of Q for in Eqs. H8I9II . as 
follows from Eq. ifnji . Matrix elements of Q are, in turn, 
the functions of eight commuting and eight anticommut- 
ing variables. Although we are interested in the highest 
order terms in anticommuting variables the calcula- 
tion of Fi [Q] and F2 [Q] is too cumbersome to be done by 
hand. The calculation can be managed most efficiently 
by employing the symbolic manipulation package epic- 
GRASS [li. The outputs of the epicGRASS (the highest 
order terms in anticommuting variables) need to be fur- 
ther integrated over all the anticommuting variables, and 
finally over all the commuting variables except 'eigenval- 
ues' [ll|. After those steps we change to the A- variables 
of Ref. jl3|, and arrive at the representation for T and 
in terms of a three-fold integral. The details of this 
procedure are outlined in the Appendix A. Here we only 
give the final expression: 



/oo /'Oo />! 

dXiJ d\2 j dAFi (A,Ai,A2) 



(ttvY 



exp {-£ (Ai A2 - A)} /i (A, Ai, A2) n (A, Ai, A2) , (17) 



/OO nOO nl 

dXi J d\2 j dAF2(A,Ai,A2) 



X exp{-e(AiA2-A)}n(A,Ai,A2)A*(A,Ai,A2), (18) 
where e = 2'KvNe, and 

1-A2 



^ (A, Ai, A2) 



(Af + A2 + A2 - 2AA1A2 - 1) 



F^ (A, Ai, A2) = 1 - A2 + (A? - l) A^ + (A^ - l) A^, 



F2 (A, Ai, A2) = 2 (1 - A? - A2 - 2AA1A2 + SA^A^)' . 

The remaining factor 11 (A, Ai, A2) contains all the infor- 
mation about the dissipation channels and comes from a 
calculation of the relevant superdeterminant (cf. for 
the GUE case) 
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(.9 + A) 



Af 



2gAiA2 + + A2 - 1 



Af 



= n(A,Ai,A2), 



A/f 



(19) 



r 



where g = (w^ + 7^) / {2'n:v^v'^) , and we have also used: 

StrQK^^Ai{Xi\2-X). (20) 

To generalize Eqs. Ijl7ll8|l to the case of non-equipotent 
dampers we just need to replace 11 (A, Ai, A2) with 



n(<?,,A,Ai,A2) = n 



i.9^ + A) 



^gf + 2.g,AiA2 + A2 + A^ - 1 



see, for example [uLEit. It can be verified that Eq. ifTzjl 
yields the same result for T as follows from adopting the 
final formula of Ref . ■ 



E. Special case of uniform damping. Comparison 
with naive calculation. 



Next, we compare results of the present (supersymmet- 
ric) calculation with the results of Ref. 0| for the case of 
uniform damping M — 0, 0. In that special case the 
naive calculation of Ref. |3] should be exact. In order to 
obtain T we need to evaluate the integral 



I{x) 



dXi 



dX exp {ix (Ai A2 — A)} 



X Fi (A,Ai, A2)/i(A, Ai,A2) 



(21) 



where we have denoted x = ie. The Fourier transforma- 
tion with respect to the x— variable, 

/•OC 

I {t) = / I (x) cxp {—ixt} dx, 



has a meaning of averaged response power in the time do- 
main for a system without dissipation. It can be written 
in a more convenient form: 



after a lengthy but straightforward procedure arrive at a 
very simple expression: 

I{t) = Aire {t) . 

which can be immediately Fourier-inverted, yielding 

This is equivalent to the first moment of the transmitted 
power given by 



(tti/) 



(23) 



and indeed coincides with the value predicted by the 
naive calculation of Ref. . 

The same steps can be repeated when calculating 
the second moment T^. One starts with Fourier- 
transforming the right-hand side of Eq. Ijl8|l . then 
changes to the variables u and z, carries out the remain- 
ing double integration explicitly and finally apphes the 
Fourier-inversion. Intermediate calculations are too long 
to be reproduced in the paper, but the final result reads: 
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— = - (5 + 28£ 



(ttz^) 



7e'] - 



(5 + 2e + e2) 



+ —El (e) (10 + lOe + Se^ + 
+ ^Ei (e)(-10+10e-3e2 + e''' 



(24) 



where 



El (z) 



-ds. 



This matches perfectly with the corresponding result of 
Ref. H. 



/ (t) = 277 



dXi 



dXn 



dXS (A-AiA2+i) 



xFi(A,Ai,A2)m(A,Ai,A2) 



(22) 



After performing A— integration, we make the change of 
variables: u = AiA2,z = Af suggested in Ref. and 



III. NUMERICAL RESULTS FOR THE 
MOMENTS OF THE TRANSMITTED POWER. 

The predictions Eqs. H17I18|I of the supersymmetric 
calculations can be compared with Monte Carlo evalua- 
tions of the first two moments of T. Towards this goal we 
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numerically generated an ensemble of N x N real sym- 
metric matrices H typically choosing 1500 ensemble real- 
izations and taking N — 1000. Theprocedure is almost 
identical to that described in Ref. The entries in H 
are constructed using a Gaussian random number gener- 
ator such that HijHki = {l/N) {SikSji + SuSjk). To sim- 
ulate the case of the uniform damping we use F = el. To 
simulate the case of a finite number of decay channels we 
take the diagonal T = diag{7, 7, 7, 0, 0} with M < 
N identical positive entries. Then, for every ensemble re- 
alization we generate the off-diagonal elements of the re- 
solvent matrix according to Gy {E) — [EI + iT — H] , 
modeling in this way the response at a site i due to ex- 
citation at the site j, with E standing for the spectral 
parameter, and both i and j chosen to be large than M 
to avoid direct coupling to the damping channels. 





Figure 1: FIG. 1. (T) and (T^) are plotted on log-scale as 
the functions of parameter e for the case of uniform damp- 
ing. The solid lines represent theoretical predictions (Eqs. 
11711 and 11811 '). For each numerically obtained ^T) and 
(r^) (represented by dots) 1500 samples of \Gij {E)\ , i ^ j 
were computed. Five sigma error bars were computed based 
on the observed variances of T and T^. 

Let us first consider the case of the uniform damping: 
r — el. For a fixed matrix size N and fixed value of 
the spectral variable E we explore a range of e. For 
E ~ the modal density dN/dE is given by v = I/tt. 
Mean level width 7 = 2itv {QEr) in this case is identical 
to e = 2TTvNe. In Fig. 1 we compare both moments 
of power T as given by Eqs. Ij2^il24ll with the results 
of Monte Carlo simulations for several values of e. It is 
evident that numerical results correspond well with the 
theoretical curves. 

To repeat the same procedure for finite number of lo- 



cal dampers M we evaluated the three-dimensional in- 
tegrals in Eqs. II17I18|I numerically for a broad range 
of the scaled mean level width 7 [1| . The difficulties of 
the numerical integration arising due to the singularity 
of /i(A, Ai,A2), are easy to overcome by employing the 
change of variables suggested in Ref. [l3|- The results 
are presented in Fig. 2 and also show a good agreement 
with the theory. 



IV. RELATIVE VARIANCE 

Dissipation within the framework of the present ap- 
proach is parametrized in terms of quantities 5, M and e. 
At the same time those quantities are not readily acces- 
sible experimentally, and in any case were not measured 
in the work reported in Ref. For this reason any di- 
rect comparison with those measurements is not feasible. 
Nevertheless, by choosing plausible values for the relevant 
parameters we can investigate the sign and magnitude of 
the discrepancies arising between the predictions for the 
relative variance of the transmitted power calculated in 
the two theories under discussion. As a result of such 
comparison we found that the discrepancy between the 
naive analysis and the present (supersymmetric) calcula- 
tion is similar to one reported previously in j3| between 
the naive theory and actual measurements. 

The comparison is carried out by first considering the 
mean spectral energy density (mean square response) 
in the time domain: E{t) = \Gij (i)|^, where Gy (t) is 
the (band-limited) time-domain Green's function. Sim- 
ilar quantities were studied previously in the context 
of the delay time distributions in chaotic cavities [Tfll |. 
Their statistics can be obtained from the inverse Fourier 
transform of the two-point correlation function T (17) = 
Gij (E) G*j [E + n) 0, 113 with respect to n. The ex- 
pression for T (fl) can be obtained by replacing 2ie with 
O in the derivation of Eq. itTzIl , (see also 113). Thus, 

E{r) - /r /r dXidX2U (r, Ai, A2) / (r, Ai, A2) 

e(AiA2-r+l)e(T-AiA2 + l)(l-(T+AiA2)^) , . 

(A5 + A^ + (T+AiA2)^-2AiA2(r+AiA2)-l)'' ' ^ ' 



where r — t/ {2t:vN) is a dimension- 
less time, /(t,Ai,A2) ^ (-^1 - l) -^2 + 
{Xl-l)Xl + 1 - (r + AiA2)', and n(r,Ai,A2) = 

(5 + 2r + A1A2)''' {g^ + 25A1A2 + A? + Ai - iy^"\ 
The generalization to the case of non-equipotent 
channels is straightforward: n(r, Ai,A2) = 

nf {g^ + 2r-f A1A2) (5? + 2g.AiA2 + A? + Al - 1)"'^'. 
The naive method yields a simpler expression for spectral 
energy density @,[3,E|: 



E{t\ 



Eo 1 



2<7T 



-M/2 



(26) 



where the initial logarithmic decay rate a is proportional 
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Figure 2: FIG. 2. Mean power (main figure) and mean square power (inset) for: (a) M — 4; (b) M — 10; (c) M = 40; (d) 
M = 400. Solid lines represent theoretical predictions (Eqs. 11711 and (O). For each numerically obtained (T) and (T^) (dots) 
1500 samples of \Gij {E)\^ , were computed. We imposed the restrictions: i 7^ j, and i > M,j > M for the non-uniform damping 
case, to avoid 'recording' the response from damped sites or from the 'source' site j, and to correspond to the assumptions in 
the theoretical analysis. For the numerically obtained mean power, twenty, ten and five sigma error bars were computed for 
the scaled mean level width 7 = 0.1, 1.0, 10.0 respectively. They were based on the observed variances of T and T^. Error bars 
for the power variance are not shown, because they are smaller than the size of the dots. For M = 4 the theoretical prediction 
for the variance does not exist. 



to the mean resonance width, given by a Porter-Thomas 
distribution. 

In Ref. i E{t) was measured experimentally, and 
fitted into the naive result Ij26|l to extract values for 
M and a. The two parameters were further used to 
predict the relative variance of T {relative variance = 
(r2) / {Tf - 1). Having the exact result (Eq. (EU) we 
can now attempt to explain the 20 — 30% over-prediction 
of relative variance reported in Ref. @|. Clearly, by spec- 
ifying certain values for M, g and e the wave scattering 
in an ergodic sample can be completely described, since 
both spectral energy density and relative variance are 
fixed uniquely. Further assuming that E{t) as given by 
Eq. II25II is the "measured" energy density of our system, 
we can repeat the procedure of Ref. . Namely, we fit 
it to E(T)naive in Order to calculate relative variance ac- 
cording to the two-parameter naive formula used in Q for 
comparison with actual measurements. Such a fit allows 
to extract values for M^aive, o'naive and Eq that may, or 
may not, correspond to the exact values. The true value 
of relative variance as determined from Eqs. Ijl7ll8|l may 



then be compared to the corresponding naive prediction. 

By a numerical three-parameter fit over the same dy- 
namic range (of a factor of e^") as in 0, we obtained 
values for Eq , anaive and Mnaive ■ In spite of the naivete 
of the model the fits were generally quite good, as ob- 
served in 0, and we can substitute the obtained values 
into the formula for the relative variance from Ref. Q 

(T^) ^ 9 M (M - 2) 

{Tf ^ ^ a {M - 4) (M - 6) 

- 4 ^1 + (T —j^^ ^2 I , (27) 

where 
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«2 



M3 

8^ 



exp 



M 
a 

M 

a 



Em 



M 
a 



trum have been calculated both analytically and numer- 
ically for ergodic dissipative structures. In the context of 
the wave scattering the former corresponds to the mean, 
and the latter to the mean square of the transmitted wave 
power T. 



Ek {z) - 




The results for several values of parameters are sum- 



M 


9 


(7 


naive 


exact 


10 


20.017 


0.497 


59.881 


59.492 


20 


20.017 


0.989 


14.419 


14.397 


20 


10.033 


1.930 


7.582 


7.574 


14 


2.918 


4.066 


4.990 


4.776 



Table I: Relative variance in absence of overall damping. 



M 


9 


e 


a 


naive 


exact 


1 


1.0 


1.0 


1.302 


7.711 


6.801 


1 


2.0 


1.0 


1.205 


7.908 


6.809 


1 


5.0 


1.0 


1.126 


8.195 


6.784 


4 


10.0 


1.0 


1.330 


7.194 


6.134 


6 


9.0 


0.5 


1.008 


10.557 


8.611 



Table II: Relative variance in presence of overall damping. 



The ergodicity assumption is implicit by virtue of our 
replacement of the actual differential operator describing 
wave motion by a large random symmetric matrix. Dis- 
sipation is taken to act both locally in space ("localized 
dampers" or dissipative channels) and uniformly within 
the sample. 

In accord with earlier results the presence of 
nonuniform, or localized, sources of dissipation requires 
the use of an elaborate nonperturbative technique — 
the so-called zero dimensional supersymmetric non-linear 
cr-model — to obtain the moments of the transmitted 
power. It is found that the naive approach fails to cor- 
rectly describe mean square power; the failure is ascrib- 
able at least in part to the assumption of real Gaussian 
eigenmodes inherent in that approach. The supersymme- 
try technique allows one to bypass the difficulty of iden- 
tifying eigenmode statistics, and to arrive at expressions 
which are in agreement with Monte Carlo simulations, 
and appear to be in better agreement with experimen- 
tally measured values of variance. 



marized in Tab. 1 and Tab. 2. It appears that in the 
absence of overall damping (Tab. 1) the actual value of 
relative variance is very close to its naive estimate. How- 
ever, when we consider the case of a small number of 
strong dampers in a system with a uniform background 
e ^ (Tab. 2), the difference becomes similar to the 
discrepancy found in A more definitive comparison 
of Eqs. (I17I18|I with measurements awaits an experiment 
in which the values of e and the gi can be ascertained 
independently. 

V. CONCLUSIONS 

In the present paper the special cases of two and four 
point correlation functions of the transmitted power spec- 
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Appendix A: EVALUATION OF THE 
SUPERINTEGRAL 

In this Appendix we elaborate on steps leading to the 
main results of Section II, expressed by Eqs. H17I18|I . We 
start by evaluating Fi [Q] and F2 [Q] with epicGRASS. 
The program extracted terms of lowest and highest or- 
der in anticommuting variables, which are, generally, the 
only terms needed. The lowest order term was found to 
be unimportant since the resulting integrands are not sin- 



gular at the boundary Then, we simplified the 

output of epicGRASS with Mathematica and reduced 
the superintegral to a multiple integral over commuting 
and anticommuting variables [ll|. 

The elements of matrix Q are introduced into the 
epicGRASS in terms of the parametrization of Ref. 

Eight commuting variables are: the "eigenval- 
ues" /ii,/i2,/i, the parameters of SU (2) group m,r,s, 
and two "angles" (^1 and (p2- The integration region in 
Eqs. Ijl5ll6|l corresponds to —00 < ^1,^2,^,7", s < 00, 

< /i < 1, < (^1, (^2 < 27r jiil. Then, after epicGRASS 
extracts the highest order term in anticommuting vari- 
ables, we have, for example, for Fi [Q]: 

Fi [Q] = — 32z^ — 32ziZ2 cos tpi cos siniy9i sin (p2 

— Zi (36 cos (fi cos (^2 + 12 cos ip^ sin ^p2 
+12 cos (p2 sin ipi+A sin ipl sin (^2) 

— z| (36 sin (fil sin + 12 cos (fi sin (^2 
-|-12 cos ip2 sin (fll + 4 cos cos 1P2) , 

where zi^2 = /ii,2^1 + Mi 2) ^-^d ^ — «/^\/l — M^- 

Eight anticommuting variables are readily integrated 
out according to the convention / dxx = 1/ (27r)^^^. 
Note, that this convention is different from the one we 
took in the beginning of Section II in the derivation of 
generating function. However, this discrepancy has no 
influence on the remaining process, as long as we use the 
integration measure of Ref . Finally, integrating over 
the "angles" as well as over the parameters of SU (2) we 
arrive at: 

Fi [Q] = -16 {zf + zl- 2z^) , (Al) 

where we indicated the integration (which does not affect 
other factors in the integrand in Eq. iflRll ) by tilde. 

Upon the substitution of "eigenvalues" into Eq. |IA1|I 
we can compare our expression for T with the final for- 
mula of Ref. • We switch to the combinations 

which are the final variables appearing in the resulting 
expression of Ref. Two results match perfectly, 

and we can proceed with the corresponding calculation 
of the second moment of the transmitted power. Before 
doing that we again change variables, this time — to the 
"eigenvalues" of Efetov's parameterization, according to 

\X^2 - A1A2 ± ^J{\i-i){xi-i). 

The domain of the integration has to be modified as well: 

1 <Ai,A2 <oo,-l <A< 1. The Efetov's "eigenvalues" 
are somewhat more convenient for the calculations done 
at the end of Section II, where we compared the exact 
and naive results for the first two moments in uniform 
damping case. 
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The analogous procedure for yields 

F2 [Q] = 4 (4x^ - Axxi - 4:XX2 + x\+ xl 
+2xiX2 + 8zl + %zl - l&z^f 



where 2:1^2 = 1 + = 1 ~ and after passing 

to Efetov's variables in Eqs. IIA1IA2II we obtain the final 
results of Section II — Eqs. II17I18|I . 

(A2) 



